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Abstract The objective of this paper is to give a rigorous analysis of a stochastic spatial model 
■ of producer-consumer systems that has been recently introduced by Kang and the author 

to understand the role of space in ecological communities in which individuals compete for 
, resources. Each point of the square lattice is occupied by an individual which is characterized 

by one of two possible types, and updates its type in continuous time at rate one. Each 
' individual being thought of as a producer and consumer of resources, the new type at each 

, update is chosen at random from a certain interaction neighborhood according to probabilities 

proportional to the ability of the neighbors to consume the resource produced by the individual 
to be updated. In addition to giving a complete qualitative picture of the phase diagram of the 
spatial model, our results indicate that the nonspatial deterministic mean-field approximation 
of the stochastic process fails to describe the behavior of the system in the presence of local 
interactions. In particular, we prove that, in the parameter region where the nonspatial model 
displays bistability, there is a dominant type that wins regardless of its initial density in the 
' spatial model, and that the inclusion of space also translates into a significant reduction of 

, the parameter region where both types coexist. 



1. Introduction 

> . 

Q>^ , To understand the role of space in ecological communities in which individuals compete for re- 

sources, Kang and the author [TU] recently introduced a stochastic spatial model of producer- 
consumer systems on the lattice. Based on their numerical simulations of this spatial model and 
their analytical results for its nonspatial deterministic mean-field approximation, they concluded 
\ that the inclusion of space drastically affects the outcome of such biological interactions. Their 

CO ' results for the models with two species are reminiscent of the ones found in [15] for the spatially ex- 

plicit Lotka-Volterra model and in |12] for a non Mendelian diploid model: in the parameter region 
where the nonspatial model displays bistability, there is a dominant type that wins regardless of 
^ ' its initial density in the spatial model, while the inclusion of space also translates into a significant 

^ ■ reduction of the parameter region where both species coexist. In the presence of three species, other 

disagreements between the spatial and nonspatial models appear. The main purpose of this paper 
is to give rigorous proofs of some of the results stated in [lOj for the stochastic spatial model with 
two species, which we now describe in detail. 

Model description. Each point of the d-dimensional square lattice is occupied by exactly one 
individual which is characterized by one of two possible types, say type 1 and type 2. Each individual 
is thought of as a producer and consumer of resources, and we described the dynamics based on 
four parameters aij > that denote the ability of an individual of type i to consume the resource 
produced by an individual of type j, that we call resource j. To avoid degenerate cases, we also 
assume that each resource can be consumed by at least one type: oij -|- a2j > 0. Each individual 
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dies at rate one, i.e., the type at each vertex is updated at rate one, with the new type being chosen 
at random from a certain neighborhood according to probabihties proportional to the abihty of the 
neighbors to consume the resource produced by the individual to be updated. More precisely, the 
state space consists of all functions rj that map the square lattice into the set of types {1,2} and 
the dynamics are described by the Markov generator 

Lfiv) = E E ,/^^^f''^L , if (v..) - fm (1) 

where fi{x,r]) denotes the number of type i neighbors 

fi{x,r]) = card{y ^ x : \yi - xi\ -\ \- \yd - Xd\ < M and 7]{y) = i} 

and where the configuration rjx^i is obtained from the configuration r/ by assigning type i to vertex 
X and leaving the type of all the other vertices unchanged. To avoid cumbersome notations, we will 
write later x ~ y to indicate that vertices x and y are neighbors, i.e., 

y^x and I2/1 - xi| + |?/2 - a;2| H V\yd-Xd\<M. 

The positive integer M is referred to as the range of the interactions. The fraction in ([T]) means 
that the conditional probability that the new type is chosen to be type i given that vertex x is of 
type j at the time of an update is equal to the overall ability of the neighbors of type i to consume 
resource j divided by the overall ability of all the neighbors to consume resource j. Note that the 
assumption on the parameters, aij + a2j > for j = 1, 2, allows to define 

ai := ail (ail + a2i)"^ G [0,1] and 02 := a22 (ai2 + a22)"^ G [0,1] 
and to rewrite the Markov generator ([1]) in the form 

indicating that our model reduces to a two-parameter system. Note also that if one of the original 
parameters in equation ([TJ is equal to zero then the denominator in 

aijfi{x,r]) _ {1 - aj) fi{x,r]) 

* 7^ J l-jj 



aij fi{x, r]) + a2j f2{x, v) «j fj{x, ??) + (1 - aj) fi{x, t]) 

may be equal to zero, in which case the numerator is equal to zero as well. When this happens, we 
have the following alternative and define ([3]) as follows: 

1. (ij = fi = 0. In this case, motivated by the fact that there is no neighbor of type i, we define 
the probability ([3]) at which vertex x becomes type i to be zero. 

2. aj = fi = 1. In this case, motivated by the fact that there are only neighbors of type i, we 
define the probability ^ at which vertex x becomes type i to be one. 

We point out that, in addition to be the most natural from a biological point of view, the two as- 
sumptions above also make the transition probability ([3]) continuous with respect to the parameters, 
which is a key to some of our proofs. 
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The mean-field approximation. As pointed out in [TO], one of tlie most interesting aspects 
of the spatial model is that it results in predictions that differ significantly from its nonspatial 
deterministic mean-field approximation, thus indicating that the spatial component plays a key 
role in these interactions. We recall that the mean- field model is obtained by assuming that the 
population is well-mixing, and refer to Durrett and Levin [8] for more details. This assumption 
results in a system of coupled differential equations for the densities of each type, and since the 
densities sum up to one, the mean-field model of ([2]) reduces to the one-dimensional system 

dui ai2 ui U2 a2i U2 ui 



dt ai2 Ui + 022 U2 ail ^^i + «21 U2 

(4) 

_ (1 - 02)^1 (1 - ui) _ (1 - ai)ui (1 - ui) 
02 (1 - ui) + (1 - 02) ui ai ui + {1 - oi) (1 - ui) 

where Ui denotes the density of type i individuals. Interestingly, the limiting behavior of the mean- 
field model depends on whether the parameters ai and 02 are smaller or larger than one half, which 
has a biological interpretation in terms of altruism and selfishness. More precisely, 

1. type i is said to be altruistic whenever Oj < 1/2 since in this case the resources it produces 
are more beneficial for the other type than for its own type, 

2. type i is said to be selfish whenever Oj > 1/2 since in this case the resources it produces are 
more beneficial for its own type than for the other type. 

The mean-field model has two trivial equilibria given by ei = 1 and 62 = 0, respectively. By setting 
the right-hand side of equation (jH equal to zero and assuming, in order to study the existence of 
nontrivial interior fixed points, that the product ui (1 — ui) 7^ 0, we find 

{l-a2){ai-l/2)ui = (l-ai)(a2-l/2)(l-ni) 

from which it follows that 

^ (l-ai)(a2-l/2) 

(l-a2)(ai-l/2) + (l-ai)(a2-l/2) 

is the unique fixed point that may differ from zero and one. To analyze the global stability of the 
fixed points, we observe that the sign of the derivative (jl]) is given by 

sign((l - a2)(ai - 1/2) ui - (1 - ai)(a2 - 1/2) (1 - ui)) 

(5) 

= sign((l - a2)(ai - 1/2) (ui - e^) + (1 - ai)(a2 - 1/2) {ui - e*)), 

which leads to the following three possible regimes for the mean-field model. 

1 - A selfish type always wins against an altruistic type. 

Proof. Without loss of generality, we may assume that ai < 1/2 < 02 indicating that type 1 
is altruistic whereas type 2 is selfish. In this case, the first line of equation ([5]) implies that the 
derivative of ui is always negative when ui G (0, 1) hence there is no fixed point in (0, 1), and the 
trivial equilibrium ei is unstable while the trivial equilibrium 62 is globally stable. □ 

2 - In the presence of two selfish types, the system is bistable. 

Proof. Since oi > 1/2 and 02 > 1/2 we have e* G (0, 1). The second line of equation ([5]) implies 
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that the derivative of ui has the same sign as ui — so the interior fixed point is unstable and 
the boundary equihbria locahy stable: the system converges to ei if the initial density of type 1 is 
strictly larger than e*, and to 62 if the initial density of type 1 is strictly smaller than e*. □ 

3 - In the presence of two altruistic types, coexistence occurs. 

Proof. Since ai < 1/2 and 02 < 1/2 the fixed point e* again belongs to (0, 1) as in the presence 
of two selfish types. However, the second line in equation ([5]) now implies that the derivative of ui 
has the same sign as e* — ui so the interior fixed point is globally stable whereas the two boundary 
equilibria ei and 62 are unstable. □ 

The stochastic spatial model. In order to state our results for the spatial model, we start by 
defining some of the possible regimes the spin system ([2]) can exhibit. 

Definition 1 - For the stochastic spatial model ([2]), we say that 

1. type i wins if starting from any configuration with infinitely many type i, 

lim P {rit{x) =i) = 1 for all x G Z^. 

t~^oo 

2. the system clusters if starting from any translation invariant configuration, 

lim P {r]t{x) = r]t{y)) = 1 for all x,y G Z'^. 

t— ^00 

3. type i survives if starting from any conGguration with infinitely many type i, 

hminf P {r]t{x) = i) > for all x G Z'^. 

4. both types coexist if they both survive. 

As pointed out above, the analytical results for the spatial and the nonspatial models differ in 
many aspects, which reveals the importance of the spatial component in the interactions. These 
disagreements are more pronounced in low dimensions and emerge when both types are selfish 
or both types are altruistic, while in the presence of one selfish type and one altruistic type the 
analytical results for both the spatial and nonspatial models agree. 

To begin with, we focus on the one-dimensional nearest neighbor process. In this case, one 
obtains a complete picture of the phase diagram based on a simple analysis of the interfaces of 
the process which almost evolve according to a system of annihilating random walks, whereas the 
more challenging study in higher dimensions and for larger ranges of interactions relies on very 
different techniques. We point out that the one-dimensional nearest neighbor process appears as a 
degenerate case in which coexistence is not possible while increasing the spatial dimension or the 
range of the interactions may result in coexistence of altruistic types. 

Theorem 2 - Assume that M = d = 1. Then, 

1. Except for the case oi = 02 = 1, the process clusters. 

2. Assuming in addition that ai < 02, type 2 wins. 
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fli = fli = 1/2 ai = 1 fli = fli = 1/2 fli = 1 

Figure 1. Phase diagrams of the nonspatial (left) and spatial (right) models. 



The exclusion of the case ai = a2 = 1 is simply due to the fact that, under this assumption, a 
vertex may flip to a new type if and only if all its neighbors already have this type, which implies 
that the process has infinitely many absorbing states. While in the absence of space one type wins 
regardless of the initial densities only in selfish-altruistic systems, Theorem[2]indicates that the most 
selfish/least altruistic type always wins in the one-dimensional nearest neighbor case. In particular, 
coexistence no longer occurs even in altruistic systems. Although this case is pathological, it is also 
symptomatic of the general behavior of the spatial model for which the existence of a dominant 
type that always wins may occur even when both types are altruistic. The next result shows that, in 
any spatial dimension and for any dispersal range, selfish-altruistic interactions result in the same 
outcome in the spatial and nonspatial mean- field models. 

Theorem 3 - Assume that oi < 1/2 < 02- Then, type 2 wins. 

Even though its rigorous proof is not straightforward, the result is predictable since, when one type 
is selfish and the other type is altruistic, regardless of the type of the vertex to be updated and 
the spatial configuration, each of the selfish neighbors is individually more likely to be selected to 
determine the new type than each of the altruistic neighbors. In contrast, in the case of selfish- 
selfish interactions covered by Theorem |4] below, the spatial and nonspatial models disagree in all 
dimensions and for all dispersal ranges. Recall that in this case the mean-field model is bistable, 
indicating that the long-term behavior is determined by both the parameters and the initial den- 
sities. The next result shows on the contrary that, including local interactions, even if both types 
are selfish, enough asymmetry in the parameters implies that the most selfish type wins. 

Theorem 4 - For all oi < 1 there is p > such that type 2 wins when 02 > 1 — p. 

Numerical simulations in two dimensions indicate more generally that, when type 2 is selfish, it wins 
whenever ai < 02- In the neutral case, the obvious symmetry of the evolution rules implies that 
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none of the types wins, but simulations indicate that the population again clusters, with boundaries 
getting sharper as the common value of ai = 02 > 1/2 increases. This suggests that, in the presence 
of two selfish types, the long-term behavior of the one-dimensional nearest neighbor process occurs 
in two dimensions as well, a property that we believe is true in any spatial dimension and for 
any range of interactions. Looking finally at altruistic-altruistic interactions, recall that coexistence 
always occurs in the mean-field approximation whereas the least altruistic type always wins in the 
one-dimensional nearest neighbor process. Our last two results show that the long-term behavior of 
the process in higher spatial dimensions and/or for larger dispersal ranges is intermediate between 
these two extreme behaviors. 

Theorem 5 - For all {M,d) 7^ (1, 1), there is p > such that coexistence occurs when 

ai < p and 02 < p- 

In words, sufficiently strong mutual altruism translates into coexistence of both types. Numerical 
simulations of the two-dimensional system suggest in addition that two altruistic types coexist in 
the neutral case oi = 02 < 1/2 therefore the coexistence region of the spatial model stretches up to 
the parameter point corresponding to voter model interactions. Interestingly, our last result shows 
that there exists a parameter region for which even altruistic types cannot coexist, thus indicating 
that the inclusion of local interactions translates into a reduction of the coexistence region in general 
and not only in the one-dimensional nearest neighbor case. 

Theorem 6 - There exists p > such that type 2 wins whenever 

ai < p and 02 > 1/2 — p. 

We note that, even though the inclusion of local interactions always translate into a reduction of 
the coexistence region, numerical simulations indicate that this reduction becomes less pronounced 
while increasing the spatial dimension or the range of the interactions. The results for the mean-field 
model and spatial model when {M,d) 7^ (1, 1) are summarized in Figured! In the phase diagram 
of the spatial model, the dashed regions are the ones that are covered by our theorems, while the 
thick lines are the critical curves suggested by simulations. For a detailed analysis of the nonspatial 
model and additional results based on numerical simulations for the spatial model in the presence 
of three species, we refer the reader to the companion paper |10j . 

The rest of this paper is devoted to the proofs. To avoid cumbersome notations and messy 
calculations, we prove all our results, with the exception of Theorem \5\ when the range of the 
interactions is equal to one. However, we point out that all the techniques used in the proofs easily 
extend to larger dispersal ranges. The reason for excluding Theorem [5] from this rule is that it is the 
only result that shows a disagreement between the one-dimensional nearest neighbor case and the 
other cases, thus indicating that the range of the interactions plays a specific role in the presence 
of strong altruistic- altruistic interactions. 

2. The one-dimensional spatial model 

This section is devoted to the proof of Theorem [2] which gives a complete picture of the long-term 
behavior of the one-dimensional process with nearest neighbor interactions. To begin with, we study 
the evolution in the non-neutral case, and look at the process starting with only type 2 individuals 
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in an interval of length + 1. The following lemma shows that the survival probability of the 
individuals of type 1 decreases exponentially with N when ai < a2- 

Lemma 7 - Assume that M = d = 1 and ai < a2- Then, 

P {r]t{x) ^ 2 for aiJ x G Z I r]o{x) = 2 for all xG[0,N]nZ) > 1 - 

where cq := (1 — 02) (1 — ai)^^ < 1. 

Proof. First, we observe that the process is attractive, which directly follows from the monotonicity 
of the probabilities in ([3]) with respect to the fractions of neighbors of each type. In particular, it 
suffices to prove the result when starting with exactly + 1 individuals of type 2. To understand 
how the number of type 2 vertices evolve, let x G Z and assume that the neighbors of x have 
different types. Then, the rate at which vertex x jumps from type 2 to type 1 is 

(1 -a2)/i(x,r/) l-a2 
" — - - 1 - 02 (6) 



(1 - 02) /i(a;,r/) + 02 /2(x,??) (1 - 02) + 02 
while the rate at which it jumps from type 1 to type 2 is 

^ (1 - ai)f2{x,v) _ 1 - Qi 

^ ■ ai /i(x,r/) + (1 - 01) /2(x,r/) ai + (1 - ai) 



1-ai. (7) 



In particular, as long as there are at least two individuals of type 2, the number of such individuals 
performs a random walk that jumps to the right at rate 2c+ and to the left at rate 2c_ since, due to 
one-dimensional nearest neighbor interactions and the choice of the initial configuration, at all times 
the set of vertices of type 2 is an interval. Moreover, on the event that type 2 survives, the leftmost 
vertex of type 2 converges almost surely to —00 while the rightmost vertex of type 2 converges 
almost surely to +00, therefore the probability to be estimated is equal to the probability that the 
asymmetric random walk described above tends to infinity. The result then reduces to a standard 
estimate for asymmetric random walks. More precisely, let Zn denote the number of individuals of 
type 2 after n updates of the system, and let 

Pi := P {Zn = 1 for some n > | Zq = i). 

Then, using a first-step analysis, for all i > 2 we have 

(c+-hc_)Pi = c+ Pj+i c_ Pi_i and 

P^+i-Pi = co{Pi-Pi-i) = ■■■ = 4-^(^2-^1) 

where cq = c_/c_|_ = (1 — 02) (1 — ai)""'^- Since Pi = 1, we deduce that Pi = Cq^""^ and 

P {rit{x) ^ 2 for all X G Z I r/o(x) = 2 for all x G [0, N] n Z) 

> P {r]t{x) ^ 2 for ah X G Z I {x : ?7o(x) = 2} = [0, iV] n Z) > 1 - Pn+i = 1 - c^- 

This completes the proof. □ 



With Lemma [7] in hands, we can now establish the first part of Theorem [2] in the non-neutral 
case as well as the second part of the theorem. 
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Lemma 8 - Assume that M = d = 1 and ai < a2. Then, type 2 wins. 

Proof. Let be a positive integer. Then, starting from any initial configuration with infinitely 
many vertices of type 2, with probability one, there exist z G Z and t < co such that 

rit{x) = 2 for x = z,z + l,...,z + N. 

Using in addition that the process is attractive, and that the evolution rules are translation invariant 
in space and time. Lemma [7| implies that 

P {rit{x) —7- 2 for all x G Z | cardja; : r]o{x) = 2} = oo) 

> P {r]t{x) ^ 2 for all X G Z I 770(2;) = 2 for aU x G [0, iV] n Z) > 1 - c^. 

Since this holds for all N and since cq < 1, the lemma follows. □ 



To complete the proof of Theorem [2l it remains to show that, when oi = 02 7^ 1, the process 
clusters. Due to the particular geometry of the configurations in the one-dimensional nearest neigh- 
bor case, the proof reduces to the analysis of an auxiliary process that we shall call the interface 
process: the spin system (^j) defined on the translated lattice D = Z + 1/2 by setting 



for all V G 



In words, the process has a one at sites f G D which are located between two vertices that have 
different types, and a zero at sites which are located between two vertices that have the same type, 
so the process keeps track of the interfaces of the spatial model. 



Lemma 9 - Assume that M = d = 1 and oi = a2 ^ 1. Then, the process clusters. 

Proof. Thinking of each site v G B as being occupied by a particle if ^(w) = 1 and as empty 
otherwise, the idea of the proof is to establish almost sure extinction of this system of particles, 
which is equivalent to clustering of the spatial model. First, the proof of Lemma [7| indicates that a 
vertex of either type, say type i, changes its type 

1. at rate one if none of its two nearest neighbors is of type i, 

2. at rate c_ = c+ > defined in ([6j) and ([7]) if its neighbors have different types. 

This induces the following dynamics for the interface process: a particle at site v G D jumps to 
each of its empty nearest neighbors at rate c+ and two particles distance one apart annihilate each 
other at rate one. More formally, the Markov generator is given by 

v& \w-v\=l 

+ E E li^(^) = ^(^)} [/(^''--) - /(^)] 

v€B \w-v\=l 

where configuration ^v(r^w is obtained from ^ by exchanging the contents of vertices v and w while 
configuration ^y^y^ is obtained by killing the particles at v and w if they exist. Since particles can 
only annihilate, the probability that a given site is occupied by a particle decreases in time so it has 
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a limit when time goes to infinity. Since in addition the evolution rules of the process are translation 
invariant, this limit does not depend on the site under consideration: 

there exists a such that lim P {^,t{v) = 1) = a for all v € B. 

t—>-oo 

In other respects, given any two particles alive at time t, at some random time s > t almost 
surely finite, either one of these particles is killed due to a collision with a third particle or both 
particles annihilate each other. The latter follows immediately from the fact that one-dimensional 
random walks are recurrent while each time two particles are distance one apart, there is a positive 
probability that they annihilate at their next jump when ai = 02 7^ 1. This implies that the limit 
a must equal zero, from which we conclude that 

lim P {7]t{x) r]t{y)) < lim P ( x + i - ^) = l) = 

t— >-oo t^-oo ^ — ^ y y Z J J 

for all x,y G Z. This completes the proof. □ 



3. Altruistic-selfish interactions 



This section is devoted to the proof of Theorem [3] which states that, in the presence of altruistic- 
selfish interactions, the selfish type always wins. The proof is divided into four steps. The first step 
is to show that, under the assumptions of the theorem, the set of type 2 dominates its counterpart 
in a process that we shall call perturbation of the voter model. In particular, in order to establish 
the theorem, it suffices to prove that type 2 wins for this new process. The reason for introducing 
a perturbation of the voter model is that, contrary to process ([2]), it can be studied using duality 
techniques, and the second step of the proof is to describe its dual process in detail, while the third 
step is to construct selected dual paths that are key to proving that type 2 survives. Finally, the 
fourth step combines these selected dual paths to a block construction in order to prove that not 
only type 2 survives but also type 1 goes extinct. Before giving the details of the proof, we note 
that the third step will be used again in the proof of Theorem [6] while the fourth step in the proof 
of both Theorems H] and El but these two steps are detailed in this section only. We also point out 
that Theorem [3] can be proved without the use of a block construction, but since it is needed in 
the proof of Theorems H] and [Gj we follow the same approach for all three theorems. 

Coupling with a voter model perturbation. We first observe that, under the assumptions of 
the theorem, there exists a constant /j > fixed from now on such that 

ai < and 02 > — ^ (8) 

in which case the set of type 2 for the spatial model (l2|) dominates stochastically its counterpart in 
a certain perturbation of the voter model that we denote later by (^t). The dynamics of this voter 
model perturbation depend on a single parameter: 

' d{l -P) + P 

and can be described as follows. As in the original spatial model, the type at each vertex x is 
updated at rate one, but the new type is now chosen according to the following rules: 
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1. with probability 1 — e > 0, the new type at vertex x is chosen uniformly at random from the 
set of the nearest neighbors, 

2. with the residual probability e > 0, the new type at vertex x is chosen to be type 1 if all the 
nearest neighbors are of type 1, and type 2 otherwise. 

More formally, the dynamics are described by the Markov generator 
= [(l-e)/i(x,e) + el{/2(x,0 = 0}] 

+ E [(i-e)/2(x,o + ei{/2(x,e)/o}] - /(e)] 

where configuration ^z,* is obtained from ^ by assigning type i to vertex x and leaving the type of 
all the other vertices unchanged. Then, we have the following result. 

Lemma 10 - Assume that r]Q = and ([8]) holds. Then, 

p {r^^{x) = 2) > P {it{x) = 2) for aiJ x G Z'^ and all t > 0. 

Proof. This follows from certain inequalities between the transition rates. When all the neighbors 
of vertex x are of the same type at the time of an update, then vertex x becomes of this type with 
probability one in both processes. Also, given that x is of type 2 and has at least one neighbor of 
each type, the rate at which it becomes 1 for process ([2|) is bounded from above by 

{l-a2)fi{x,ri) ^ {1- p)fi{x,'n) 



(l-a2)/i(x,r/)+a2/2(x,??) (1 - p) r/) + (1 + p) /2(x, r/) ^^^^ 

[1- p) + 2pf2{x,r]) d{l-p)+p 

which is the rate at which it becomes 1 for process (jlOp . Similarly, given that x is of type 1 and 
has at least one neighbor of each type, the rate at which it becomes 2 for process ([2]) is 

(1 - ai)/2(x,ry) ^ {1 + p) f2{x,ri) 



ai fi{x,'n) + (1 - ai) /2(x,r/) (1 - p) fi{x,r]) + (! + />) /2(x,??) 

(12) 

(l + p)/2(x,??) d{l- p) f2{x,r]) + p 

which is the rate at which it becomes 2 for process (jlOp. To see that the last inequality is indeed 
true, we introduce the functions 

[1 - p) + 2p z d{l - p) + p 

and observe that we have the equalities 

1 + P 



^(1) = h{l) = 1 and 9{{2dy^) = h{{2d)-^) 



2d{l - p) + 2p 



Since in addition g is concave on (0, 1) whereas h is linear, the last inequality in (I12p follows. The 
inequalities on the transition rates in (jlip and (jl2p give the lemma. □ 
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Duality with branching coalescing random walks. In view of Lemma \T0\ standard coupling 
arguments - see Liggett [13j, Section II.l, for details about coupling techniques - imply that, to 
prove the theorem, it suffices to prove that type 2 wins for the process (jlOp . Recall that the reason 
for introducing this process is that, contrary to the original spatial model, it can be studied using 
duality techniques. First, we construct the process ([10]) graphically as follows: 

1. for each oriented edge e = {x ,y) ^ 1/ x with x ~ y 

(a) we let A(x,y) be a Poisson point process with intensity {2d)~^{l — e), 

(b) we draw an arrow from vertex y to vertex x at times t G A(j;, y) to indicate that the 
individual at x mimics the individual at y. 

2. for each vertex x G Z'^ 

(a) we let A(x) be a Poisson point process with intensity e, 

(b) we draw a set of 2d arrows from each y ~ x to vertex x at times t G A(x) to indicate 
that X becomes of type 1 if all its neighbors are of type 1, and of type 2 otherwise. 

We refer the reader to the left-hand side of Figure [2] for an illustration of the graphical representa- 
tion. The type at any space-time point can be deduced from the graphical representation and the 
configuration of the system at earlier times by going backwards in time. Declare that a dual path 
from space-time point (x, T) down to (y, T — s) exists, which we write (x, T) ^ (y,T — s), whenever 
there are sequences of times and vertices 

so = T-s < si < ••• < Sm+i = T and xq = y, xi, . . . , x^ = x 

such that the following two conditions hold: 

1. for i = 1, 2, . . . , m, there is an arrow from Xj_i to Xj at time Sj and 

2. for i = 0, 1, . . . , m, there is no arrow that points at the segments {xj} x (sj, Sj+i). 
The dynamics imply that vertex x is of type 1 at time T if and only if 

^o(y) = 1 for ah y G Z'^ such that (x,r) i (y,0). (13) 

Note that the set-valued process 

6(x, T) = {y G Z'^ : (x, T) i (y, T - s)} for ah < s < T 

defines a system of branching coalescing random walks where particles jump at rate 1 — e to one of 
their nearest neighbors chosen uniformly at random and split at rate e into 2d particles which are 
sent to each of their neighbors. In addition, whenever two particles are located on the same vertex 
they coalesce. See the left-hand side of Figure [2] for an illustration where the system of random 
walks is represented in thick lines. This process is reminiscent, though not identical, of the dual 
process of the biased voter model [1] and we also refer to [U [11] for similar dual processes. 

Selected dual paths and random walk estimates. To bound the probability that a given 
space-time point (x,T) is of type 1, we now construct a dual path {Xt{x)) that keeps track of a 
specific particle in the dual process. For each i = 1,2, . . . ,d, we define the hyperplane 



Hi := {z = {zi,Z2,...,Zd) eR'^ : z, = 0} 
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as well as the deterministic times Tj = iciN where ci = 4e~ . Recall that the constant e is defined 
in equation Q above. Also, we let dist(3;, H) denote the Euclidean distance between a vertex x and 
its orthogonal projection on a set H. The dual path starts at Xq{x) = x and is defined recursively 
from the graphical representation as follows. For alH > 0, define 

s{t) = mi{s > t :T - s £ A{Xt{x)) ov T - s £ A{Xt{x),y) for some y ~ Xt{x)}. 

For all s G {t,s{t)), we set Xs{x) = Xt{x) while at time s{t) we have the alternative: 

1. If T — s{t) £ A{Xt{x), y) for some y ~ Xt{x) then Xg^i-^{x) = y. 

2. IfT - s{t) G A(Xt(x)) and t G {Ti^i,Ti) for some i = 1, 2, . . . , d, then 



well-defined. In words, the dual path travels backwards in time down the graphical representation 
following the arrows from tip to tail. When a set of 2d arrows is encountered, the process jumps 
in the direction that makes it closer to the hyperplane Hi between time Tj-i and time Tj while 
after the last time Td the process jumps uniformly at random in all directions. This, together with 
the properties of the graphical representation introduced above, implies that the selected dual path 
always jumps at rate one. At each jump, with probability 1 — e, the target site is chosen uniformly at 
random from the set of the nearest neighbors whereas, with probability e, the process moves so as 
to decrease its distance to Hi until time Ti when a similar mechanism operates in the direction of 
the second axis, and so on. In order to later compare the process pop with oriented site percolation, 
we start by collecting important properties of the selected dual path. These properties are based 
on random walk estimates and are given in the following three lemmas. 

Lemma 11 - There exist Ci < oo and 71 > such that, for all i = 1,2, . . . ,d, 



3. IfT 



Xsit){x) Xt{x) and dist{Xg(^t){x), Hi) < dist{Xt{x), Hi). 
s{t) £ A{Xt{x)) and t > then 




P{dist{XT,{x),Hi) > N/2 I dist(XT,_i(x),F,) < AN) < Ci exp(-7iiV) 



for all N sufficiently large. 




(l-6)+6 



lim h-^E{dist{Xt+h{x),Hi)-dist{Xt{x),Hi)) 



= —e. 



(14) 



In particular, recalling the definition of Tj and Tj_i, we obtain 



E{dist{XT,{x),H^) I dist{XT,_,ix),H^) = m) 



m-eiTi-Ti.i) = 0. 



The lemma then follows from large deviation estimates for the Poisson distribution and the Binomial 
distribution together with standard coupling arguments. □ 
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Lemma 12 - For any constant C > Ade'^ and all i = 1,2, . . . ,d, 

P {dist{Xt{x),Hi) > N for some t E (T^, CN) \ dist{XT,{x), Hi) < N/2) < 3 exp(-\/iV) 

for all N sufRciently large. 

Proof. Denote by {Z^) the discrete-time one-dimensional simple symmetric random walk and 
introduce the discrete random variable Jg that counts the number of times the ith coordinate of 
the process {Xt{x)) jumps by time s. Since the ith coordinate evolves according to the continuous- 
time simple symmetric random walk run at rate d~^{l — e) after time Tj, we have 

P{dist{Xt{x),Hi) > N for some t G {Ti,CN) \ dist{XT,{x), Hi) < N/2) 

< P {JcN > No) + P {Zn i {-N,N) for some n < No \ Zo € {-N/2, N/2)) 

for all A'^o > 0. Taking No = d^^CN, large deviation estimates for the Poisson distribution imply 
that the first term on the right-hand side can be bounded by 

PiJcN>No) < PiJcN > d-\l-e)CN) < C2 exp{-j2N) < exp{-VN) 

for suitable constants C2 < 00 and 72 > 0, and all N sufficiently large. To estimate the second 
term, we use the reflection principle and Chebyshev's inequality. For all 9 > 0, we have 

P{Zn i {-N,N) for some n < A^o I ^0 G (-iV/2, iV/2)) 

< PiZn^ {-N/2, N/2) for some n < iVo I ^0 = 0) 

< 2 P(ZjVo i i-N/2,N/2) \ Zo = 0) < 2 exp{-9N/2) E{exp{9ZNo) \ Zo = 0) 

< 2 exp{-9N/2) [E{exp{9Zi) \ Zq = 0)]^" < 2 exp{-9N/2 + Noln(l){9)) 

where the function <j3{9) is the moment generating function of Zi. Since 

= ln(^^^^^) = In (^1 + y + 0(^2)^ < 9^ 

when 9 > is small enough, taking 9 = 4/\/iV, we can conclude that 

P{Zn i {-N,N) for some n < A^o I ^0 G (-iV/2, iV/2)) 

< 2 exp(-6'iV/2 + iVoln</'(0)) < 2 exp(-2^/iV + 16d~iC) < 2 exp(-ViV) 

for all N sufficiently large. This completes the proof. □ 

Lemma 13 - Let C3 > 4de~^ and x G (-2iV, 2A^]"'. Then, for all N large, 

P {Xt{x) i (-4iV, ANf for some t < C3N or Xc;n{x) (/ {-N, Nf) < 7d exp(-\/iV). 

Proof. Since dist{Xt{x), Hi) is stochastically smaller than the rate one simple symmetric random 
walk on the set of nonnegative integers with a reflecting boundary at zero and starting from the 
same initial state, the proof of Lemma [12] directly implies that 

P {dist{Xt{x),H^) > AN for some t < C3N) < 3 exp(- ViiV) (15) 
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for all N large and alH = 1, 2, . . . , d. Moreover, by Lemmas 1111 and 1121 

P{dist{Xc3N{x),Hi) > N I dist{Xt{x),H,) < AN for all t < C3N) 

< P (dist(XT, (x), H,) > N/2 I dist(XT,_i ix),Hi) <m) 

(16) 

+ P{dist{Xt{x),Hi) > N for some t £ (T„ CgiV) | dist{XT,ix), Hi) < N/2) 

< Ci exp(-7iiV) + 3 exp(-^/iV) 

for all N large and alH = 1, 2, . . . , d. Finally, combining (jlSp and (jl6p . we conclude that 

^ (-4iV,47V]'^ for some t < C3N or Xcsn{x) (/ (-iV,iV]^) 

< P{Xt{x) ^ (-4iV,47V]'^ for some t < C3N) 

+ P{Xc:,n{x) i {-N,NY and Xt(x) G (-47V,4iV]'^ for all t < C3N) 

< 3dexp{-VAN) + dCi exp(-7iiV) + 3d exp{-VN) < 7d exp{-VN) 
for all N sufficiently large. This completes the proof. □ 

Block construction. To complete the proof of Theorem [3l the last step is to use a rescaling 
argument, a technique also known as block construction, which has been introduced by Bramson 
and Durrett [3]. The general idea of the block construction is to couple certain good events related 
to the interacting particle system properly rescaled in space and time with the set of open sites of 
an oriented site percolation process on the lattice 

7i = {(z, n) = (zi, . . . , n) G Z"' X Z : zi + . . . + + n is even and n > 0}. 

We refer to Section 4 in Durrett [7\ for more details about this technique and a rigorous definition 
of oriented site percolation. Lemmas II 111131 give us the appropriate space and time scales to compare 
the process (llOp with oriented percolation. More precisely, we let T = C3N where the constant C3 
is as in Lemma [T3l and declare {z,n) G ?^ to be good whenever the event 

n{z, n) = {Ut{x) = 2 for ah x£Nz+ (-iV, Nf n Z^} 

occurs. Then, we have the following lemma. 

Lemma 14 - For all N sufRciently large, 

P{Xt(x) i {-AN,mY for some {x,t) G {-2N,2NY x (0,r)) 

+ P ((±ej, 1) is not good for some i = 1, 2, . . . , d | (0, 0) is good) 
< 7d{{AN)'^ + 1) exp{-^/N) 
where Ci denotes the ith d-dimensional unit vector. 

Proof. By construction, there exists a dual path (x,T) | {Xt{x),0). Therefore, Lemma [T3l and 
the duality relationship ()13p imply that 

P ((±ej, 1) is not good for some i = 1,2, . . . , d \ (0, 0) is good) 

= P {£,t{x) = 1 for some x G (-2iV, 2N]'^ \ ^o{y) = 2 for all y G (-A^, Nf) 

< P {^o{Xt{x)) = 1 for some x G {-2N, 2N]'^ \ ^o{y) = 2 for all y G {-N, iV]^) 

< (47V)'^ sup,e(-27V,2iV]<* P iUXrix)) = 1 | ^y) = 2 for all y G {-N,N]^) 

< (AN)'' ^^P.e(~2N,2N]' P{XT{x)i{-N,NY) < 7d(47V)^exp(-^/]V) 
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for all N large. Lemma [13] also implies that for all x G {—2N, 2N]'^ we have 

P{Xt{x) i (-4iV,4Af]°' for some t G (0,r)) < Id exp(-ViV) 

for all N sufficiently large. The lemma follows. □ 

The right-hand side of the inequality in the statement of Lemma [14] tends to zero, therefore both 
probabilities on the left-hand side tend to zero, as tends to infinity. Convergence to zero of the 
second probability together with the translation invariance of the evolution rules of the process in 
space and time imply that, for all 7 > 0, there exists large such that 

P ((z ± ej, n -|- 1) is good for alH = 1, 2, . . . , d | (z, n) is good) > 1 — 7. (17) 

Convergence to zero of the first probability in the statement of the lemma implies that, with 
probability close to one when N is large, the good event in (|17p only depends on the realization 
of the graphical representation in a finite space-time box which, together with (|17p . implies the 
existence of events Gz^n such that 

1. Gz^n is measurable with respect to the graphical representation in 

{Nz (-47V, MSlf} X [nT, (n l)r] for ah (z, n) e 

2. For all 7 > there exists large such that P {Gz,n) > 1 — 7, 

3. We have the inclusions of events 

Gz^n n ri(2;,n) C ri(2; lb ej, n -|- 1) for all i = 1, 2, . . . , 

These three conditions are the assumptions of Theorem 4.3 in Durrett [7], from which we deduce 
the existence of a coupling such that the set of good sites contains the set of wet sites of a certain 
oriented site percolation process with parameter 1 — 7 and finite range of dependency when the 
scale parameter is large. Since such a percolation process is supercritical when 7 is small enough, 
starting with infinitely many vertices of type 2, we have 

liminf P (^i(x) = 2) > for ah x G 'L'^ . 

This only proves survival of the type 2 individuals, and the last step of the proof is to also establish 
extinction of the type 1 individuals, which is not obvious since oriented site percolation has a 
positive density of closed sites. The existence of an in-all-directions expanding region void of type 1 
is proved in the following lemma which we will apply again to prove Theorems [1] and [UJ 

Lemma 15 - For all x € Z,'^ we have P (Ct(x) = 1) — t- as time t — t- 00. 

Proof. The proof follows an idea of Durrett [6j which relies on the lack of percolation of the dry 
sites, where a site in the percolation process is said to be dry if it is not wet. First, we define the 
oriented graphs Gi and G2 with common vertex set Ti and respective edge sets 



El = {{{zi,ni),{z2,n2)) € H X Ti : zi - Z2 € {±ei, . . . ,±ed} and rii = n2 - 1} 

E2 = El U {((zi, ni), (z2, "-2)) & H X H : zi - Z2 e {±2ei, . . . , ±26^} and ni = 71-2}. 
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Given a realization of the percolation process in which sites are open with probability 1 — 7, we say 
that there is a regular dry path from (0, 0) to {z, n) if there is a sequence 

(zo,?T-o) = (0,0), (2:1,711), {zk,nk) = {z,n) £n 

such that the following two conditions hold: 

1. {{zi,ni), {zi+i,ni+i)) £ El for aU i = 0, 1, . . . , A; - 1 and 

2. the site (zj, n^) is dry for all i = 0, 1, . . . , /c. 

Extending condition 1 by replacing Ei with the larger edge set E2, we call the sequence a generalized 
dry path. Assuming that initially only site (0, 0) is closed, we introduce 

1. Dq{Qi) = sites that can be reached from the origin by a regular dry path and 

2. Dq{Q2) = sites that can be reached from the origin by a generalized dry path. 

We call these sets the regular dry cluster and the generalized dry cluster, respectively. Returning 
to the particle system, in order to prove extinction of type 1, we first modify our definition of good 
sites: a site (z, n) £ T-l is now said to be good if the event 

n{z, n) = {^^(x) = 2 for ah s £ [nT, (n + 1)T] and aU x£ Nz+ (-iV, N]'^ n Z"'} 

occurs. The same arguments as in the proofs of Lemmas [TTIfT^ imply that 

P ((±ei, 1) is not good for some i = 1, 2, . . . , d | (0, 0) is good) < C4 N'^ exp(-\/]V) 

for some suitable C4 < 00, thus there again exists a coupling such that the set of good sites 
contains the set of wet sites of an oriented site percolation process with parameter 1 — 7 provided 
the parameter N is large. The reason for introducing this new definition of good sites is that the 
presence of a type 1 individual at vertex x G Nz + {—N, N]'^ at time nT now implies the existence of 
a generalized dry path from some {zo,0) to {z,n). In particular, the problem of extinction reduces 
to proving the lack of percolation of the generalized dry sites when 7 > is small enough. To do 
this, the first step is to prove the existence of C5 < 00 and 75 > such that 

P {Do{Q2) t ([-m, mf X M) n 7^) < C5 exp(-75m) for 7 > small. (18) 

This is essentially Theorem 5 in Durrett [6] which states the result for the regular dry cluster but 
his proof easily extends to our case. To begin with, we define the collection of cubes 

Q = {r G : Ei=i,2,...,d+i kil < 1} and Q{z,n) = (z,n)+Q for (z, n) G ^ 

and orient the faces of the cubes by assigning the value +1 to the d + 1 faces at the top of each 
cube, and the value —1 to the d + 1 faces at the bottom. Then, we define the contour of the regular 
and generalized dry clusters as the algebraic sum of all the faces of the cubes that make up 

'^oiG^) = U {Qiz, n) : (z, n) G Do(a^)} for i = 1,2. (19) 

To prove (jlSp . the first ingredient is to observe that each contour has an equal number of plus and 
minus faces and sites below the plus faces must be closed. This holds regardless of the geometry 
of the underlying oriented graph so Lemma 4 in [6J directly applies: if the contour has length m 
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then there is a set of at least m{2{d + 1))~^ sites that must be closed. The second ingredient is 
given by Lemmas 5-6 in [B] which state that there exist Cg < oo and /i > such that the number 
of contours with m faces is bounded by /i"*. The result is proved for the regular dry cluster but 
again it is straightforward to extend the proof to the generalized dry cluster. To see this, the key is 
simply to observe that any cube that makes up the dry region ([T9|) associated with the generalized 
dry cluster has at least one corner in common with another such cube. Then, following the lines of 
the proof of Theorem 5 in [6] and using the previous two estimates, we obtain 

oo 
k=m, 

from which (|18p follows for 7 < To deduce from (jlSp the lack of percolation of the dry 

sites, we denote by Wn and Wn the set of wet sites at level n in the coupled percolation processes 
having the same sets of open sites except at level where 

Wo = {0} and Wq = {z e Z'^ : zi + Z2 -\ h is even}. 

The proof of Lemma 8 in [6J gives the existence of a constant a > such that 

Wn n ^2(0, an) C Wn for all n large on the event {Wn / for all n > 0}. (20) 

Following the proof of Lemma 11 in [6j, and using (|18p and (|20p . we obtain 

limm->-oo P (a generalized dry path starting at some {z, 0) 

intersects (-62(0, na/2) x {n}) n Ti for some n > in) = 

on the event {Wn 7^ for all n > 0}. Through the coupling with the process (llOp . we deduce the 
existence of an in-all-directions expanding region containing only good sites, thus the existence of 
an in-all-directions expanding region which is void of type 1. □ 

From Lemma [T5l we deduce that type 2 wins for the process (|10p therefore the same holds for 
the process ([2]) according to Lemma \T0\ which completes the proof of Theorem [3l 



4. Selfish-selfish interactions 



This section is devoted to the proof of Theorem |4] which claims the existence of cases of selfish- 
selfish interactions for which the most selfish type outcompetes the other type whenever starting 
with infinitely many representatives. The key to our proof is to first show that type 2 wins for the 
process with parameters ai < 1 and 02 = 1, that we denote by (f/t), and then use a perturbation 
argument to extend the result to a neighborhood of this parameter point. First, we note that, under 
the assumption 02 = 1, the dynamics of the process ^ reduce to 

Lflfifi) = ^ l{/2(x,ry)=0} -/(??)] 

I (1 -ai)/2(x,f?) (21) 

aifi{x,7]) + {l-ai)f2{x,7]) 

Invasion of type 2 individuals in the particular case when 02 = 1 is established in the following 
lemma through a comparison with the Richardson model introduced in 
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Lemma 16 - Let 5 > 0. Then, for all N sufRciently large, 

P ivcrNix) = 1 for some x e (-2iV, 2N]'^ \ fio{x) = 2 for all x G {-N, N]'^) < 5 
for some constant Cy < oo that does not depend on 6. 

Proof. First, we observe that, given that vertex x has at least one neighbor of type 2 at the time 
it is updated, it flips from type 1 to type 2 with probabihty 

(1 - ai)f2{x,fi) ^ ^ _^ 1 -Qi ^ Q 

ai fi{x,ff) + {I - ai) f2{x,ff) ~ ' fli (2(i - 1) + (1 - ai) 

while it flips from type 2 to type 1 with probability zero. This implies that, for every connected 
component A C Z'^ containing at least two vertices, if all vertices in A are of type 2 at some time, 
then this remains true at any later time. In particular, the set of type 2 individuals dominates 
stochastically the set of occupied sites in the Richardson model {^t), i-e., the contact process with 
no death, with birth parameter e{2d)~^ > under certain initial conditions. More precisely, both 
processes can be constructed from the same graphical representation in such a way that 

P i^t C {x : fjtix) = 2} for alH > Uo C {x : f?o(x) = 2}) = 1 

whenever the initial set of type 2 individuals is connected and has at least two vertices. This, 
together with the Shape Theorem for the Richardson model [16J, implies that there exists a finite 
constant C7 < oo that only depends on the parameter e such that, for all 6 > 0, 

P {fjc7N{x) = 1 for some x G {-2N, IN]'^ \ fio{x) = 2 for all x E {-N, N]'^) 
< P{ic,N 75 (-2iV,2iV]'^ I Co = {-N,N]'') < S 

for all sufRciently large. This completes the proof. □ 

Extinction of type 1 individuals stated in Theorem [4] follows from Lemmas [15] and [16] through 
the comparison of the process properly rescaled in space and time with a certain oriented site per- 
colation process in which dry sites do not percolate. Let T = C^N where C7 < 00 is the constant 
fixed in Lemma [TBI and declare {z,n) E to be a good site whenever the event 

n{z, n) = {fis{x) = 2 for ah s E [nT, (n + l)r] and ah x£Nz + {-N, N]'^ D Z°'} 

occurs. In the particular case when ai < 02 = 1, individuals of type 2 with at least one neighbor of 
type 2 cannot change, therefore it follows from Lemma [TBI that, for all S > 0, 

P {{z ±ei,n + 1) is good for alH = 1, 2, . . . , d | (z, n) is good) > 1 — 5 

for all large. Theorem 4.3 in Durrett [7] again gives the existence of a coupling such that the set 
of good sites contains the set of wet sites of a certain oriented site percolation process with finite 
range of dependency and parameter 1 — 5. To complete the proof of Theorem [H we show extinction 
of individuals of type 1 when 02 < 1 is close enough to 1 relying on a perturbation argument. 

Lemma 17 - Let oi < 1. Then, there exists p > small such that 

lim P (??j(x) = 1) = for all x E Z*^ whenever 02 > 1 — p. 
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Proof. First, we fix 7 > such that (jlSh in Lemma [TSl holds, and set 6 = 7/2. Then, we apply 
Lemma [TU] to obtain the existence of a scaling parameter N such that 

P{fis{x) = 1 for some (x, s) S {-2N,2NY x [T,2T] \ 

f/,(x) = 2 for all G (-iV,iV]'^ X [0,r]) < 5. ^^^^ 

The scaling parameter being fixed, since the transition rates are continuous with respect to the 
parameters, there exist a small p = p{N) > and a coupling of ([2]) and (I2ip such that 

P{ris{x) / f]s{x) for some {x,s) G {-2N,2N]'^ x [T,2T] \ 

rj^{x) = f]s{x) = 2 for aU {x,s) G {-N,N]'^ x [0,T]) < 5 ^^^^ 

whenever a2 > 1 — p- By combining (j22p and (|23p . we obtain 

P{r}s{x) = 1 for some (x, s) G {-2N,2NY x [T, 2r] | 

= 2 for ah (x,s) G {-N^N^ x [0,r]) < 25 = 7 

whenever 02 > 1 — p. Hence, a new application of Theorem 4.3 in Durrett gives the existence of 
a coupling such that the set of good sites contains the set of wet sites of an oriented site percolation 
process with parameter 1 — 7 when 02 > 1 — p. Having such a coupling, extinction of type 1 
individuals directly follows as previously from the proof of Lemma [T5j □ 



5. Strong altruistic-altruistic interactions 



This section is devoted to the proof of Theorem [5] which indicates that, in contrast with the one- 
dimensional nearest neighbor model, coexistence occurs in any other spatial dimensions and for any 
other range of interactions provided cooperation is strong enough. The first step of the proof is to 
show that, excluding the one-dimensional nearest neighbor case, the critical value of the threshold 
contact process is strictly smaller than one, and then use a standard coupling argument to deduce 
the existence of a small constant /) > such that type 1 individuals survive for the process with 
parameters oi = and 02 = p- The coexistence result will then follow from symmetry and the 
monotonicity of the survival probabilities with respect to the parameters. Recall that the threshold 
contact process with parameter a is the spin system [^t) in which each occupied vertex becomes 
empty at rate one, and each empty vertex with at least one occupied neighbor becomes occupied 
at rate a. Formally, the dynamics are described by the Markov generator 

hf{i) = «i{/i(^'e)/o} + E [/(e-.o) - /(e)] 

where empty and occupied vertices are identified to states zero and one, respectively. The reason 
for looking at the threshold contact process is that, when the birth rate a = 1, it is closely related 
to the process ([2]) with parameters oi = 02 = 0. 

Lemma 18 - Assume that {M,d) / (1,1)- Then, 

Uc = inf {q : P (^t / for aii t > | Co = {0}) > 0} < L 
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Proof. The rate at which a vertex becomes empty is constant while the rate at which a vertex 
becomes occupied is a nondecreasing function of the configuration, where the partial order on the set 
of configurations is the inclusion. This indicates that the threshold contact process is an attractive 
spin system. In particular, Theorem 2.4 in Bezuidenhout and Gray [2] implies that 

{a : P{^t^ for aU t > | = {0}) > 0} is an open interval. (24) 

In addition. Lemma 2.1 in [9j, which relies on a result of Liggett [14], implies that 

P (^t / for ah t > Uo = {0}) > for (M, d) / (1, 1) and a > 1. (25) 

Combining (f24l) and (|25]) gives 1 G (ac, oo) therefore ac < 1. □ 

Relying on Lemma [T8l we now prove the existence of p > small such that type 1 individuals 
survive for the process with parameters ai = and 02 = p, that we denote by (%). Note that, for 
this choice of the parameters, the dynamics reduce to 

L,f{fi) = E J, — [fifi.,i)-fm 

X&1 (26) 

+ i{h{x,f))^Q] [fifj,,2)-fm- 

Lemma 19 - For all (M, d) ^ (1,1), there is p > such that 

liminf P {fjt{x) = 1) > for all x G Z'^. 

Proof. Thinking of vertices of type 1 as occupied by a particle and those of type 2 as empty, the 
expression of (p6|) indicates that particles die at rate at most one, while empty vertices with at 
least one occupied neighbor become occupied at rate at least 

/ N ^ 1 - P ^ 1 - P 

"^^^ l-p + p((2M + l)'^-l) l + ((2M + l)'^-2)p 

provided p > is small enough since Oc < 1 according to Lemma [181 Therefore, for p > small 
such that the previous inequality holds, the set of type 1 vertices dominates stochastically the set 
of occupied sites of a supercritical threshold contact process. Hence, it suffices to prove the result 
for the latter: starting with infinitely many occupied vertices, 

liminf P {^t{x) = 1) > for ah x e Z'^ and all 

This is a direct consequence of the complete convergence theorem for the threshold contact pro- 
cess established in Cox and Durrett [5j which implies that, starting with infinitely many occupied 
vertices, the process converges in distribution to the limit obtained when starting with all vertices 
occupied, which is a nontrivial translation invariant measure when a > ac- □ 

To deduce Theorem [5] from Lemma [T9l we first invoke the monotonicity of the survival proba- 
bilities with respect to the parameters to obtain that, for p > as in the lemma, type 1 survives 
when the parameters are ai = 02 = p. By symmetry, individuals of type 2 survive as well, and so 
both types coexist, for these parameters. Finally, invoking again the monotonicity of the survival 
probabilities, we deduce that coexistence occurs for all 01,02 < p. 
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Figure 2. Pictures related to the proof of Theorems [3] and |6l 



6. Strong altruistic-weak altruistic interactions 

This section is devoted to the proof of Theorem [6] which gives the existence of a parameter region in 
which coexistence occurs for the nonspatial model but not the spatial model. The proof combines 
ideas from Sections [3] and [H First, we study the process with ai = and 02 = 1/2 using techniques 
with are reminiscent of duality techniques, though the process does not have a dual process at this 
parameter point. This is done through the new concept of breaking point that we introduce in this 
paper. Using breaking points, we construct selected dual paths as in Section [3l and using a block 
construction, we deduce that type 2 wins when oi = and 02 = 1/2. To complete the proof of the 
theorem, we rely as in Section |4] on a perturbation argument. 

Breaking points. First, we assume that ai = and 02 = 1/2. Under this assumption, the type 
at each vertex x is updated at rate one according to the following rules. 

1. If vertex x is of type 2 at the time of the update, then the new type is chosen uniformly at 
random from the set of the 2d nearest neighbors. 

2. If vertex x is of type 1 at the time of the update, then it switches to type 2 whenever at least 
one of its nearest neighbors is of type 2. 

More formally, the dynamics of the process with ai = and 02 = 1/2, that we denote by (f/t), is 
described by the Markov generator 

Lfififi) = E + E l{/2(x,??)/0} [/fe,2)-/(r?)]. (27) 

We construct the process ([27|) graphically as follows: 

1. for each oriented edge e = (x, y) G Z*^ x Z"^ with x ~ y 
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(a) we let A(x,y) be a Poisson point process with intensity {2d)~^, 

(b) we draw at all times t G A(x, y) a solid arrow y — t- x and an additional 2d — 1 dashed 
arrows from the other neighbors of x to vertex x to indicate that, if x is of type 2 then 
it mimics vertex y whereas if it is of type 1 then it switches to type 2 whenever at least 
one of its neighbors is of type 2. 

Definition 20 ~ We call {x, t) where t G A(2;, y) a breaking point if 

sup{A(y,x) n (0,t)} = sup {(A(x, z) UA(y,z)) n (0,t) : z ~ X or z ~ y} 
and the set A(y, x) n (0, t) is nonempty. 

In words, (x,t), t G A(x,y), is a breaking point if the next solid arrow pointing at either vertex 
X or vertex y we see going down the graphical representation is directed from x to y. We refer to 
the right-hand side of Figure [2] for an example of realization of the graphical representation along 
with the breaking points which are represented with black dots. The reason for looking at breaking 
points is that they allow, as in the proof of Theorem [3l to construct selected dual paths with a drift 
towards regions void of type 1 individuals and thus to prove invasion of type 2. This good property 
of the breaking points is expressed in the following lemma. 

Lemma 21 - Assume that (x, t) is a breaking point. Then 

fit{x) = 1 if and only if fit{z) = 1 for all z ~ x. 
Proof. Assume that t G A{x,y) and let s = sup{A(y,x) n (0,t)}. Then, 

f]s{x) = 2 =^ ^s{y) = 2 because there is a solid arrow x — )• y at time s 

=^ Vr{y) = 2 for all r G since no arrow points at {y} x {s,t) (28) 

=^ f]t{x) = 2 because there is a solid arrow y — )• x at time t. 
From (f28]l . we deduce that 

f]t{x) = 1 =^ f]t{x) = 1 and f]s{x) = 1 

=^ f]t{x) = 1 and fit~{x) = 1 because no arrow points at {x} x 
=^ fit{z) = 1 for all z ~ X 

where the last implication follows from the evolution rules of the process. This proves the necessary 
condition. The sufficient condition is trivial and we point out that it is always true provided an 
arrow points at (x,t) whereas the reverse holds because (x,t) is a breaking point. □ 

Motivated by Lemma [2T| we now construct from the graphical representation a certain system 
of branching coalescing random walks that we shall call a dual process as it allows to obtain infor- 
mation about the type of a given vertex, say x, at the current time, say T, based on the configuration 
at earlier times. To define the dual process starting at (x,T), the first step is to remove from the 
graphical representation all the dashed arrows that point at a space-time point which is not a 
breaking point. Then, we say that there is a dual path from {x,T) down to (y,T — s), which we 
again write (x,T) I {y,T — s), whenever there are sequences of times and vertices 

so = T-s < si < ■■■ < Sm+i=T and xq = y, xi, . . . , x^, = x 



such that the following two conditions hold: 
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1. for i = 1,2, . . . ,m, there is an arrow from to Xi at time Si and 

2. for i = 0,1, . . . ,m, there is no arrow that points at the segments {xj} x (sj, Sj+i). 

We point out that condition 1 above must hold after removal of the dashed arrows that are not 
directed to a breaking point. The dual process starting at (x, T) is the set-valued process 

fi,{x, T) = {y G : {x, T) i{y,T - s)} for all < s < T. 

See the right-hand side of Figure [2] where the dual process is represented in thick lines. One of 
the keys to proving Theorem [6] is the following duality relationship which is simply obtained by 
applying Lemma [21] at each of the breaking points that occurs along the dual process: 

Vt{x) = 1 only if f]T~s{y) = 1 for all y G r)s(x,T). (29) 

Selected dual paths. To bound the probability that space-time point (x, T) is of type 1, we follow 
the same approach as in Section [3] and construct a process {Yt{x)) that keeps track of a specific 
particle in the dual process. Recall that 

Hi := {z = {zi,Z2,... ,Zd) GM"^ : Zi = 0} for i = 1, 2, . . . , d. 

In addition, we let Tj := i C2N for i = 0,1, . . . ,d, where 

C2 := 4e-^ with e := 2d (1 - e-^/^)^ e'^^'^"^) > 

and where is a scaling parameter that will be fixed later. The choice of e will become clear in 
the proof of Lemma [22] below. The dual path starts at Yo{x) = x. To construct this path at later 
times assuming that it is defined up to time t, we introduce 

s{t) = mf{s>t:T-se A{Yt{x),y) for some y ~ Yt{x)}. (30) 

For all s £ {t,s{t)), we set Ys{x) = Yt{x) while at time s{t) we have the alternative: 

1. If (yj(x), T — s{t)) is a breaking point and t £ {Ti^i, Ti) for some i = 1,2, . . . ,d, then 

Y,^t){x) ^ Yt{x) and dist{Y,^t){x), Hi) < dist{Yt{x), Hi). 

2. Otherwise, ^^^((^(x) = y is the neighbor for which the infimum in (|30p is reached. 

In words, the process goes backwards in time down the graphical representation following the solid 
arrows from tip to tail except when a breaking point is encountered in which case the process 
crosses the arrow, either dashed or solid, that makes it closer to the hyperplane Hi between the 
times Tj„i and Tj. After time Td only the solid arrows are crossed. The process is thus embedded 
in the dual process which, together with the relationship (j29p . implies that 

fjj,(x) = 1 only if Yrix) = 1. (31) 

Motivated by condition pip , the next step is to prove the analog of Lemma [13] for this new selected 
dual path. More precisely, we have the following lemma. 

Lemma 22 - Let x G (-2A^, 2N]'^. Then, there exists Cg < 00 such that 

P {Yt{x) i (-4A^,4iV]^ for some t < CsN or Yc,n{oo) ^ {-N,Nf) < Cg exp(-79\/iV) 
for suitable Cq < 00 and 79 > 0, and all N sufficiently large. 
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Proof. To begin with, we call s G (0, T — 1) a breaking time whenever 

1. there is a unique neighbor y of vertex Ys{x) such that there is a solid arrow y — ?• Ys{x) in the 
time interval {T — s,T — s — 1/2), 

2. there is a unique solid arrow Ys{x) — )■ y in the time interval {T — s — 1/2, T — s — 1), 

3. no other solid arrows point at Ys{x) or y in the time interval (T — s,T — s — 1). 

Note that these three conditions imply that the dual path crosses exactly one breaking point 
between dual times s and s + 1. Moreover, a simple calculation shows that 

P (s is a breaking time) = 2d{l- e'^^'^f f^r'^d-'^) = g > o. 

Since in addition disjoint parts of the graphical representations are independent, it follows that the 
discrete-time process (y„(x) : n G (0,T) n Z) is a Markov chain that satisfies 

£;(dist(y„+i(x),/7i)-dist(yn(x),F,)) 

< e £' (dist(y„+i(x), //j) — dist(y„(x), ffj) I n is a breaking time) < — e 

for all n G {Ti-i,Ti) Z. This is the analog of p4p . In particular, the rest of the proof follows from 
the exact same arguments as in the proofs of Lemmas [TT]|131 □ 

Block construction. To complete the proof, we let T = CgN and declare site {z, n) £ T-L tohe a 
good site whenever the following event occurs: 

n{z,n) = {??s(2;) = 2 for all s G [nT, (n + l)r] and all X G iVz+ (-iV,iV]'^nZ'^}. 

Lemma [22] and (I3ip imply that, for all 5 > 0, 

P {{z ± ei,n + 1) is good for alH = 1, 2, . . . , d | (z, n) is good) > 1 — 5 

provided the scale parameter is sufficiently large. This, together with the exact same arguments 
as in the proofs of Lemmas 1141 and I15| implies that type 2 wins for the process (j27p . Finally, using 
the continuity of the transition rates with respect to the parameters and the same approach as in 
the proof of Lemma [TTl we deduce the existence of p > small such that type 2 wins as well for 
the process ([2]) when ai < p and a2 > 1/2 — p. This completes the proof of Theorem [6l 
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